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I. INTRODUCTION 

The knowledge of correlations is important and in- 
teresting for any system. Looking from the practical 
point of view correlation means additional information: 
if two quantities are correlated, the knowledge of one of 
them implies certain information about the other one. 
In physical systems correlations usually indicate inter- 
actions between parts of the system. The prototypical 
example is given by the Ising spin system where the 
nearest-neighbors interactions induce long-range correla- 
tions leading to a phase transition. 

The situation in random graphs is somewhat different. 
It is known that for random geometries even in the ab- 
sence of any explicit terms inducing interactions between 
vertices their degree may be correlated [H-Il]- Moreover, 
those correlations are long-range, i.e., they fall off as 
some power of distance They are generated by 

model constraints rather than by explicit interactions. It 
should be also stressed that the distance dependent cor- 
relation functions in the ensemble of random graphs are 
much more complicated objects than their fixed lattice 
counterparts [H, [1] . To see that let us take some generic 
correlation function on random graphs 



(1) 



where qi denotes the degree of the vertex i, i.e., the num- 
ber of branches emerging from it. A and B are some 
arbitrary functions depending on the vertex degree and 
|i — j'l is the graph (geodesic) distance between vertices i 
and j . For random geometry it makes no sense in general 
to choose two fixed points — that is why we sum over all 
the pairs of points on the graph. The graph distance is 
the length of the shortest path between those two vertices 
and as such it is dependent on the whole graph. That 
means that the above expression is not a two-point func- 
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tion but a highly nonlocal object. That is a fundamental 
difference between random and fixed geometries. 

In this paper, we study in detail correlations between 
degrees of vertices as a function of distance. We consider 
an ensemble of all labeled trees T(y) with a fixed number 
of vertices V, on which we define the probability measure 



(2) 



ieT 



riy denotes the partition function of this ensemble (nor- 
malization factor) and qi is the degree of vertex i; Wq's 
{q > 0) are some non- negative numbers (weights). This 
is a maximal-entropy ensemble with a given degree dis- 
tribution (see Appendix). An important property of the 
measure ([2]) is that it factorizes into a product of one- 
point measures, so it does not introduce any explicit 
correlations. This means that any observed correlations 
arise from the fact that we consider a specific set of 
graphs and not from the measure itself. 

We show that the connected degree-degree correlations 
are not zero and fall off with the distance as 



7(1) = TTgAl) - ^,(1)^1) 

(g-2)(r-2) 



[2 + {{q^)-A){l-l)] 



(3) 



Here 7fg^r(0 the joint probability that two vertices dis- 
tance I apart will have degrees q and r, respectively. 
Those correlations are disassortative. The average de- 
gree of the distance / neighbors of a vertex with degree q 
decreases 



k{q) 



{q')-4 



q + m-m-iy 



(4) 
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For I — 1 this reduces to the results obtained in Ref. . 
In the following sections we provide the detailed defini- 
tions of the quantities introduced above and derive those 
results. We will also discuss what happens for the scale- 
free trees when (g^) diverges. 

The paper is organized as follows: Sec. HIl introduces 
some basic definitions concerning correlations in random 
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trees. Then we derive the vertex degree distribution us- 
ing the field theory approach in Sec. Illll and proceed on in 
Sec. IIII Al calculating the distance dependent correlation 
functions. Two examples of Erdos-Renyi and scale-free 
trees are given in Sees. IIIIBI and IIII C[ respectively. In 
the following Sec. IIII Dl the results for scale- free trees are 
verified using Monte Carlo (MC) simulations. Final dis- 
cussion and summary of our results are given in Sec. IIVI 



Finally we define average degree of the vertices at the 
distance / from a vertex of degree q as follows: 



where 
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(13) 
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II. CORRELATIONS 



For each graph we introduce 



ijeG 
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which is the number of pairs of points with degrees q 
and r separated by the distance I. We define two further 
quantities: the number of pairs at the distance I with one 
end point of specified degree 



(6) 



and the number of all pairs of vertices at the distance I 
n(0^5]n,.,,(0. (7) 



If we want to define the joint probability TTq^riJ) we have 
two obvious choices. The first one is 



T^q.ril) 



n{l) 



(8) 



where the subscript denotes that we restrict the average 
to the ensemble of graphs for which n{l) is not zero. The 
second possibility is to use 



- n\ — V^q,r 
^q,r(l') = 



(0) 



nil)) 
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In Ref. Q we have argued that the first quenched defi- 
nition is more natural in the context of random graphs. 
However, it is much more difficult to work with. In this 
paper we will assume that the ensemble of generic trees is 
self-averaging and the two above definitions are equiva- 
lent. For a more detailed discussion of this issue we refer 
to [l]. Similarly we define 



_ K(0) 



(10) 
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and the connected two point probability 

^™"(Z) = 7f,,,(/)-^,(07f,(Z). (11) 

We further define the connected correlation function 



7r|°;^(/)^E<Zr^— (0. 



(12) 



III. GENERIC RANDOM TREES 

We consider an ensemble of all labeled trees with the 
probability measure The partition function fly is 
defined as the sum of the weights of all the trees in the 
ensemble 
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T£T{V) ieT 
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The partition function of the corresponding grand- 
canonical ensemble is defined by the discrete Laplace 
transform 



w = E 



(16) 



v=i 



We will use the field theory approach to calculate it [14 
We define the function 



W{^i) 



d(j) exp 



N 



E 



Wa 



(17) 

Its formal perturbative expansion in e will generate 
Feynman's diagrams with desired weights and symme- 
try factors (for an introduction see any textbook on field 
theory, e.g., Refs. [3 and [ll| or Ref. [13). This ex- 
pansion will, however, contain all the graphs including 
those which are not connected or contain loops. We can 
restrict the expansion to connected graphs only by con- 
sidering the function log W{fi). To obtain just the tree 
graphs we will use the expansion in N~^. According to 
Feyman's rules for the expression ([T7|) each edge in the 
graph introduces a factor and each vertex a factor 
N which together contribute N~^^^ , where E is the 
number of edges in the graph. If L is the number of in- 
dependent loops in the graph then E — V = L — \,so the 
first term of the expansion will group graphs with 
no loops, the second one graphs with one loop, and so 
on. 

That means that the contribution of tree graphs is 
given by the first term in the saddle-point approxima- 
tion. The saddle-point equation is 
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FIG. 1: Graphical representation of Eq. (|19p . Each gray bub- 
ble corresponds to the sum over planted trees given by the 
partition function Z{ii). Branches without a vertex (small 
empty circle) at one end denote a stem. 



We will denote by Z{fi) the solution of the above equation 
and rewrite it as 




where 



Z(m) ' 



(20) 



FIG. 3: Graphical representation of the partition function 
0,{q, fJ,) given by Eq. (|27|l for q = 4. Each gray bubble cor- 
responds to the sum over planted trees given by the parti- 
tion function Z{fi); the smaller black circle represents the 
root which contributes a weight factor e~'^W4; the additional 
factor comes from the fact that the relative position of 
branches in the compound tree is irrelevant. 



does not have any solution for /i < /xq. The function 
Z(/L() has a singularity at /zq given by the condition for 
the minimum 



Inserting Eq. ([T^ into Eq. ([T7)) and taking the logarithm 
to keep only connected graphs we obtain 



2F(Zo) = ZoF'(Zo), and = log 



(24) 



oo ^ 

ri(^) = e^^ ^ —jZ^{fji) Z'^(p)- (21) this singularity the partition function behaves like 



It is easy to check that 



Z{^) K, Zq- Zx\J\i - /io + ■^2(m - Mo), 



(25) 
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Figure [T] shows a graphical interpretation of Z{p): it is 
the partition function of the ensemble of planted trees 
[23} . Planted trees are the trees with a stem attached to 
one of the vertices. Its properties and resulting critical 
behavior were calculated in Refs. (T3l - [l5| . The model 
has two main phases. In the so called generic or tree 
phase the the function F[Z)/Z'^ has a minimum inside 
its domain (see Fig. [2]) and so Eq. (IT9)) . which can be 
rewritten as 



regardless of the form of the weights Wq . In this contri- 
bution we will limit our self to this phase only. Inserting 
the expansion ((25|) into Eq. (fT9|) and expanding to the 
order A/i = /i — /zq {Z2 cancels in the resulting equation) 
we obtain 



^ F"{Z,)Zl 
2 F(Zo) 



1. 



(26) 



F{Z) 
Z^ ' 



(23) 




The vertex degree distribution of this model was cal- 
culated using the correspondence with the balls in boxes 
model in Ref. [l6j . Here we rederive it using a differ- 
ent method which can be easily extended to the case of 
two-point correlations studied in Ref. [3]. Let us denote 
by fl{q; fi) the partition function of the rooted grand- 
canonical ensemble of trees with the condition that the 
degree of the root is q. Then 



' dw{q) 



--^z^(m). 



(27) 



The graphical interpretation of this equation is shown in 
Fig. [3] The sought degree distribution is proportional 
to the canonical partition function D,v{q)- Inserting the 
expansion (|25|) into Eq. (|27|) we obtain 
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FIG. 2: Graphical representation of Eq. 
larity of Z{ii). 
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FIG. 4: Graphical representation of the partition function 
Qi{q,r;fj,) given by Eq. p3|l for q — 4 and r = 3. Gray 
bubbles correspond to the partition function Z{^) and the 
smaller black circles mark the vertices with degrees q and 
r; double circles represent the I — 1 vertices along the path 
connecting them in which we sum over all possible insertions 
of the Z{^) function (see Fig. [5J. 



The last expression has a known inverse Laplace trans- 
form 

^-.v^ ^ 1 « (29) 

so finally keeping only the first terms in the expan- 
sion and fixing the normalization we obtain the formula 



7r(q) 



1 Wc^Zl 



FiZo) {q-iy. 



(30) 



Using the above formula we can give an interpretation of 
the right-hand side of Eq. 



2 FiZo) 2 F{Zo 



- 1 



9=1 



(31) 



Here we have used the fact that on trees the average 
degree {q) = 2 (in the large V limit). Please note that 

(9') -4 =((9- 2)2) >0. (32) 
The (g^) is equal to 4 only in the wi — > limit. 

A. Correlations in generic random trees 

We proceed as in the previous section but this time 
we introduce a partition function Q.i{q,r]^) of all the 
trees with two points marked, such that the points are 
at the distance I and have degrees q and r, respectively. 
Because we are considering the trees there is exactly one 
path linking the two marked vertices (see Fig. U). As in 



FIG. 5: Double circles denote the summation over all possible 
insertions of the Z{fi) function (gray bubbles). 



the previous section we can express the partition function 
nl{q,r■^i) by Z(^) I, [T^ 



^i{q,r\lj) 



(9-1)! 

k=2 



(r- 1)! 



(fc-2) 



(33) 



The last term comes from the vertices along the path for 
which we have to sum up all the possible insertions of 
the Z{ijl) function (see Fig. [5]). This summation can be 
done in the following way: 



,t^(9-2)! 



Z'i- 



dZ^^iq- 1) 

-M d F{Z) 
dZ Z ■ 



(34) 



Differentiating Eq. (|19p with respect to fi we come to the 
relation 



\ dZ Z ) 



Using it we finally obtain 



E' 



^ (g- 2)1 

q=2 ' 



z'i- 



Z'{p)- 



(35) 



(36) 



Inserting into Eq. (j33p first the above formula and then 
the expansion (I^Sj) we get 



^i{p,q\y) 



Z\ 



{q-\)\{r-\)\ 



q+r-2 



.^0 



l-2^VAAi) .(37) 



This can be further approximated by 



^i{v,q\p) 



(q-l)! (r-1)! " 

(9+r-2)+2|a(;_i)[yS7I 



X e L^o 



, (38) 
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and using Eq. ([29)1 we obtain to the leading order in V 



{Uqr) (X ^l{p,q; V) 



q+r-2 



(q^iy.ir-iy. 

Zjq Zji 



Finally, we get 

7f,(/) ^Tlq 



(g + r-2) + ((g2)-4) - 1) 
2+{{q^)-A){l-l) 
q+{{q^)-A) {l-l) 



2+((g2)_4) [l-l) 



(39) 

(40) 
(41) 



Inserting this into Eqs. pT|) and ([TH)) we obtain the re- 
sults and Q. Summing up Eq. Q over q and r we 
get the connected correlation function 



T^q.f (0 



[2 + ((g2)-4)(;-l)]^ 



(42) 



B. Example 1 

In the first example we put Wq — 1, so all the trees in 
the ensemble have the same weight. In this case F{Z) = 
Ze^. The solution of Eq. ([M]) is = 1 from which 
follows: 



1 1 



1 q^ 



'~ e{q-l)\ 
leading to 



q—1 ^ ' 



(44) 



l(<Z-2)(r-2) 1 1 



e2 (? + l)2 (g- 1)! (r- 1)!' ^^^^ 



and 



C. Example 2: Scale- free trees 



(46) 



In this example we choose Wq = q~^{q — 1)! which 
corresponds to the planar graphs studied in Refs. [l^ and 
[l5l |. Then F{Z) is given by the polylogarithm function 
Up{Z) 



(47) 



q=l 



and Eq. ([Ml) takes the form 2Li^(Zo) = IAp_i{Zq). It 
has the solution for jS < Pc with Pc ~ 2.4788 given by 



2C{fic) = CiPc - !)• At the critical value of ^ = the 
partition function no longer scales as in Eq. and in 
principle we cannot use the Laplace transform Eq. (|29|) 
anymore. However, as shown in Ref. [ISj the large V 
behavior is not changed and we expect our formula to 
hold in the large V limit. From Eq. ([30l) we read-off the 
degree distribution 



7r(g) = 



F{Zo) 



(48) 



At the critical value of /3, = 1 and the vertex degree 
distribution is scale free. The average 



{q^)^2 + 



1 



[Li^-2(^o)- Li/3-1 (^0)], (49) 



Li;3(2'o) 

diverges as /3 — /3c- Formula (jS)) leads for ^ > 1 to 
hm ^,=7(0=0. 



(50) 



This would imply that the correlations vanish in the large 
V limit. However, this limit (I50p is not uniform. It is easy 
to check that the integrated correlation functions do not 
disappear 



hm n^^d) 



1 



(l-l) 



and 



lim ki{q) = 2+——. 

V— >-oo t — i 



(51) 



(52) 



Please note that the above results are universal and valid 
for any kind of scale-free trees with /3 < 3. 

D. Monte Carlo simulations 

The results obtained in the previous sections are valid 
only in the strict V ^ 00 limit and it is clear that for 
finite V our formulas will not hold for any /. Defining 
the average distance on a graph 



(0 



^ n(0 

2^ y2 



(53) 



we may expect that the formulas will be valid only for 
Z <C (0- The scaling of (1) with the graph size depends 
on the Hausdorff's dimension dn 



(54) 



For generic trees considered here dn = 1/2. For scale- 
free trees considered in example 2 we expect 



7 



7 



f3c 



/3c -1' 



(55) 



which gives rf/f « 3. In the case of scale- free trees the 
volume dependence manifests itself by the cutoff in the 
degree distribution tt, as well [T9| . 
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FIG. 6: Average degree of distance I neighbors ki{q) for Z = 2 
to 6 and trees with 64 000 (empty symbols) and 128 000 (filled 
symbols) vertices. Each symbol denotes different l\ straight 
lines are the predictions given by Eq. (|52[). 



FIG. 7; Connected correlation function 7r|°"(/) for scale- free 
trees with 64 000 (empty circles) and 128 000 (filled diamonds) 
vertices. The solid line represents the prediction given by 
Eq. dsn. 



Expecting finite size effects to be more severe in the 
scale-free trees, we checked the V dependence perform- 
ing MC simulations of the ensemble described in the ex- 
ample 2. We have used an algorithm similar to "baby- 
universe surgery" [20| . The basic move consisted of pick- 
ing up an edge at random and cutting it. Then the 
smaller of the two resulting trees was grafted on some 
vertex of the bigger one. The most time consuming part 
of the algorithm was to find which tree was smaller. To 
save time the two trees were traversed simultaneously un- 
til one of them was filled completely. Additionally, to pick 
the attachment point from the bigger tree efficiently, the 
vertices of the trees were marked during the traversal. 
This move was supplemented with moves consisting of 
cutting up leaf nodes and attaching them to some other 
parts of the tree. This was much faster as it did not re- 
quire traversing the tree. However, the autocorrelation 
time for such moves alone was much higher, especially 
for the scale-free trees. Because those trees are at the 
phase transition between the generic and the crumpled 
phase [13 the autocorrelation time is high even for the 
tree grafting algorithm. 

We have simulated trees of the size up to 128 000 ver- 
tices. To verify to which extent the ensemble is self- 
averaging we measured the quenched quantities which 
we then compared to our predictions. In Fig. [5] we have 
plotted the measured 

as a function of q for various values of I. Please note 
the large finite size effects for Z = 2. This is to be ex- 
pected: For finite V the (g^) is also finite and actually 
grows slowly with V [loj . For larger I the agreement with 



our results is quite good. Figure [7] shows the quenched 
correlation function 

which is also well reproduced by our results. 

IV. SUMMARY AND DISCUSSION 

The appearance of the long-range correlations in 
generic trees is puzzling. Usually we expect the pow- 
erlike (scale-free) behavior to be manifested in systems 
at the criticality. The trees studied here apart from the 
scale-free example are, however, not critical. The free- 
energy density can be calculated in the infinite volume 
limit and remains an analytic function of the weights Wq 
(isl [l6t . It has been also shown that the critical behav- 
ior in random trees is not associated with the diverging 
correlation length [2l|. The correlations described here 
are thus of the structural and not of dynamical origin. A 
possible mechanism explaining it was proposed in Refs. 

and Q: in connected graphs vertices of degree one 
must have neighbors of degree greater than one. It re- 
mains, however, to be understood how this effect can be 
propagated to larger distances. 
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Appendix: Maximal entropy 

For each choice of the weights Wq and a given number 
of vertices V the ensemble ^ has a weh defined degree 
distribution Tr{q). Asymptoticahy for large V this distri- 
bution is given by Eq. (jSOj. Ref. f25] contains a proof 
that the probability measure ([2]) has maximal entropy 
among all the measures producing the distribution T:{q). 
Here we repeat their arguments for completeness. 

We start with an expression for the entropy plus the 
necessary Lagrange multipliers to force the constraints 

'5 = E PiT)\ogP{T) + A ^ p(T) 



Ter 



9=1 ' \Ter 



{q)V . (A.l) 



Differentiating the above with respect to P{T) we get 

logP(T)-A-l + EVg(r), (A.2) 



leading to 



P(T) 



.A,",(T) 



(A.3) 



Putting 



— , and e^" 



Wo 



(A.4) 



we obtain Eq. We will now prove that this measure is 
a unique solution of Eq. (|A.ip satisfying the constraints, 
at least in the y — > oo limit. Let us assume that we 
have another set of weights Wq that produces the same 
probability distribution Tr{q) ([501 



T:{q) = 



1 



WqZl 



1 



WqZo 



F{Zo){q-iy. F{Zo){q-l)y 
It follows that: 



F{Z^) fZo\ 

Wa = TTTTTT- 1^1 Wq, 



F{Zo) \Zo 



(A.5) 



(A.' 



hence 



n-. = n 



FjZo) 
F{Zo) 



Wa 



FjZo) 
F{Zo) 



Zo 



n^,.-(A.7) 



But this gives identical probability measure to Eq. ^ 
because of the condition X^ieT 9^ ~ — 2 valid for each 
tree T. 
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